Anomalous Hall Effect in Ferromagnetic Semiconductors in the Hopping Transport 

Regime 
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We present a theory of the Anomalous HaU Effect (AHE) in ferromagnetic (Ga,Mn)As in the 
regime when conduction is due to phonon-assisted hopping of holes between localized states in the 
impurity band. We show that the microscopic origin of the anomalous Hall conductivity in this 
system can be attributed to a phase that a hole gains when hopping around closed-loop paths in 
the presence of spin-orbit interactions and background magnetization of the localized Mn moments. 
Mapping the problem to a random resistor network, we derive an analytic expression for the macro- 
scopic anomalous Hall conductivity . We show that is proportional to the first derivative 
of the density of states Q[e) and thus can be expected to change sign as a function of impurity band 
filling. We also show that depends on temperature as the longitudinal conductivity axx within 
logarithmic accuracy. 

PACS numbers: 75.50.Pp,73.50. Jt,72.20.Ee 

Diluted magnetic semiconductors (DMS) are very in- 
teresting new materials, still not understood in full de- 
tail. Although it is well accepted that the origin of ferro- 
magnetism in Gai-^Mna; As, the most thoroughly studied 
DMS, is a hole- mediated effective interaction between the 
localized S = 5/2 Mn spins, the appropriate description 
of this interaction, and a closely related question of the 
nature of hole carriers, are still controversial issues. Sam- 
ples with the highest Curie temperatures show metal- 
lic behavior j3] and are well described by the valence- 
band hole-fluid model 0- On the other hand, for low 
Mn concentrations, i.e. when x is less than about 0.03, 
(Ga,Mn)As is insulating and the hole transport likely oc- 
curs in the Mn impurity band j^. The low-temperature 
conduction in this case is by phonon-assisted hopping. 

Anomalous Hall Effect (AHE) has been an important 
characterization tool for itinerant ferromagnets |^ 
and played an important role in the experimental study 
of DMS ferromagnetism as well • Experimentally AHE 
is manifested as an additional term in the Hall resistivity 
of the sample, usually assumed to be proportional to the 
magnetization: 



when conduction is by hopping between strongly local- 
ized states in the impurity band, rather than through 
extended valence-band states. Our approach is based 
on a modified Holstein's theory of the Hall effect in hop- 
ping conduction (for related work on AHE in colossal 
magnetoresistance manganites see 0|). 

We start from a model of spin-3/2 hole carriers, hop- 
ping between randomly distributed localizing centers 
(Mn acceptors) in the mean spin-splitting field of 5 = 5/2 
Mn local moments. The corresponding Hamiltonian is: 



Eh 

ia(3 



(1) 



g^y = RoB + RsM. 



Standard theories of AHE in metallic ferromagnets 
attribute it to the modification of impurity scattering 
in the presence of spin-orbit interactions. Recently, how- 
ever, it was realized that a purely geometrical mechanism 
of AHE is possible Hi, [lOlllll. In particular, a k-space 
Berry phase theory of AHE [g has been proposed, and 
applied specifically to (Ga,Mn)As in the metallic trans- 
port regime. This theory is very successful in describing 
AHE in metallic (Ga,Mn)As, and we believe that the 
Berry-phase mechanism of AHE should be relevant in 
general for all itinerant ferromagnets. 

In this Letter we propose an analogous geometrical- 
phase theory of AHE in (Ga,Mn)As, but in the case 



(2) 

Here i,j label localized hole states at Mn acceptor sites, 
a,P — 3/2, . . . , — 3/2 are spin-3/2 indices, are ran- 
dom on-site potentials, which model the random poten- 
tials localized holes feel due to Coulomb interactions with 
nearby charged Mn acceptors and As antisite defects, h 
is the spin-splitting mean field due the p-d exchange in- 
teraction between the hole- and localized Mn-spins, t is 
the spin-3/2 operator and Uajis is a spin-dependent hop- 
ping amplitude which we assume to be of the following 
form: 



t 



iaj(3 



(3) 



where ^ is the localization length and Aij — (e/2c)B • 
{vi X Tj) is a phase factor due to the external magnetic 
field B (we use h = 1 units). In most of the following 
we will omit this factor, t = diag{t^/2i'ti/2iti/2i'td,/2) is 
a diagonal hopping matrix in the case when fij is along 
the spin quantization direction and TZij rotates the spin 
quantization axis from z to fij . The hopping parameters 
t3/2 and ti/2 can be calculated microscopically 0. The 
fact that 7^ ^1/2 is a direct consequence of the spin- 
orbit interactions. 

We assume that the Fermi energy lies in the region 
where all states are strongly localized and thus the only 
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possible conduction mechanism is through the phonon- 
assisted hopping. Thus hole-phonon interaction must be 
exphcitly included in the model: 



q 

The full Hamiltonian is then given by 

H = Hh + Hh-ph + Hph- 



(4) 



(5) 



Local transport properties of the system described by 
Eq. Q can be most conveniently calculated using a com- 
bination of the usual linear-response theory and a per- 
turbation expansion in powers of tiajp/\ei — ej|, which is 
justified by our assumption that all states are localized 
[3 Holstein (IJ was first to realize that to cap- 

ture the (ordinary) Hall effect under hopping conduction 
conditions, one needs to consider processes where ampli- 
tudes of direct and indirect (through a third site) hops 
between two sites interfere. Such processes appear at 
third order in the perturbation theory in tiaj0/\ei — ej|, 
while the longitudinal hopping conductivity can be de- 
scribed by second-order terms. To capture AHE, one 
needs to consider third-order terms as well. Evaluating 
the linear response expression for the current lij between 
two localized states [l5|, including terms of up to third 
order in the hopping amplitudes and up to second order 
in hole-phonon interaction, one obtains: 



GijVij + ^ Fijk{Vik + Vjk)- 

k 



(6) 



Here Vij —Vi — Vj is the electrochemical potential differ- 
ence between the sites i and j. Gij is an effective conduc- 
tance between a pair of sites, which at low temperatures, 
such that \ei — ej\/T ^ 1, is given by: 



a/3 



i. .„|2g-e»c,73/T 

^lQ.Jp\ C- 



(7) 



Here 70 is a constant, proportional to the phonon density 
of states and the square of the hole-phonon interaction 
matrix element, and we have introduced the following 
notation: 



1 



(8) 



We have assumed in the above that h = hz. Henceforth 
we will also assume that h is sufhciently small, so that we 
can restrict ourselves to linear response effects with re- 
spect to h. In particular this means that we disregard the 
influence of spin splitting on the longitudinal conduction. 
We expect this to be a reasonable approximation at low 



Mn concentrations, when the magnetization is relatively 
small. 

The first term in Eq. © thus represents the usual map- 
ping of the hopping conduction problem onto a random 
resistor network 0, 0, 0| . The second term in Eq. ((BJ 
contains the physics responsible for AHE. The quantity 
Fijk is given, in the low temperature limit, by the follow- 
ing expression: 



F, 



8T ^ 



^(\<^h-,\-^io,k-y-<^h-,j[3) /T 



(9) 



Let us discuss the physical meaning of this equation. The 
process, that leads to the appearance of the term, pro- 
portional to Fijk, in the expectation value of the current 
between sites i and j, is a quantum interference process, 
where probability amplitudes of a direct hop between 
any two sites in a triad (i, j, k) and an indirect hop (i.e. 
through the third site) interfere. Such a process requires 
the participation of at least two phonons, hence the factor 
7q. The crucial component of Fijk is the imaginary part 
of the product of three complex spin-dependent hopping 
amplitudes taken along the sides of the triangle formed 
by the sites i, j, k. It is directly related to the phase that 
a hole gains while hopping along a closed-loop three-site 
path. It is this phase that ultimately leads to the ap- 
pearance of the macroscopic anomalous Hall conductiv- 
ity. Note that if the mean spin splitting field h is taken 
to zero, Yl^ai^^^iUaj/Btj/sk-rtk-yia) vanishes identically for 
any triad {i,j,k). It is the presence of a macroscopic 
magnetization in the sample, or, in other words, broken 
time reversal symmetry, that results in a nonzero contri- 
bution from this term. Another important property of 
the third-order term in the current expectation value is 
that it leads to the appearance of a current between the 
sites i and j, that is driven not by the potential differ- 
ence on the bond {ij), but by the potential differences on 
bonds {ik) and {jk). This property implies that the cur- 
rent, proportional to Fijk, is on average transverse to the 
direction of the local electric field, as one would expect 
for the Hall current. 

Eq. © reveals the microscopic origin of the AHE in the 
system. To obtain the anomalous Hall conductivity of a 
macroscopic sample, one needs to find a certain char- 
acteristic microscopic triad that determines the macro- 
scopic conductivity. Since the distributions of both dj 
and Fijk are very broad due to their exponential depen- 
dence on the random on-site potentials and Mn acceptor 
positions, the characteristic triad can not be obtained 
by a simple averaging over all possible triads. Instead, 
one needs to find a relationship between this generalized 
hopping conduction problem and a certain percolation 
problem, exactly like in the case of usual hopping con- 
duction [17I llij. It is, however, not an easy task, due 
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FIG. 1: A backbone node with an associated triad. / is the 
current flowing in a macrobond terminating at site 1 in re- 
sponse to the external electric field E. 



to a tensor nature of the transport brought about by the 
second term in Eq.©. We believe that the best approach 
to this problem is the one used in Refs. 19, 20, ,21] to in- 
vestigate the problem of ordinary Hall effect in hopping 
conduction regime. This approach starts from the as- 
sumption that the problem can be solved perturbatively 
in the second, Hall term in Eq.(|HJ. This is certainly a 
valid assumption if one wants to find the Hall conductiv- 
ity to first order in the magnetic field (or magnetization 
in our case). Thus one imagines the following two-step 
iterative procedure. Initially, the second term in (jSJ is 
put to zero and external electric field is applied to the 
network. One then solves Kirchhoff's equations 



(10) 



and finds voltages at all nodes. The external electric 
field is then turned off and the zeroth-order voltages are 
substituted in the Hall term of Eq.®. The Kirchhoff's 
equations at the second iteration take the following form: 



as well). The longitudinal conductivity of the network 
is then given by axx ~ Gc/L, where L is the correlation 
length of the cluster. The correlation length depends on 
temperature only as a power law |22| and we will not 
calculate it exphcitly. Following Refs. 0,113,13, we 
now assert that the Hall conductivity in the system is as- 
sociated exclusively with the current carrying backbone 
of the percolating cluster. The backbone can be visual- 
ized as a supernetwork of macrobonds and nodes, with 
typically three macrobonds intersecting at a node. The 
average distance between the nodes is L. Macrobonds 
can be assumed to be simply one-dimensional chains of 
conductances [l^]. It is then clear that the exponen- 
tially dominant Hall current sources will be located at 
the nodes of the percolating cluster backbone. Since the 
nodes are far away from each other, we may regard them 
as independent sources of the Hall emf. Thus the prob- 
lem reduces to finding the Hall emf developed at a single 
isolated node of the backbone in response to the cur- 
rent flowing through the node [ij, |2fl|, |2l| . Applying the 
above described two-step iterative procedure to a single 
node with an associated triad of sites shown in Fig^ we 
obtain: 



V 



H 



32 



G13G32 + G12G23 + G31G12 
^-^"123 



(13) 



G13G32 + G12G23 + G31G12 

where is the Hall emf, / is the current flowing through 
the node, and we have neglected unimportant numerical 
prefactors. Since the average distance between the nodes 
is L, the macroscopic Hall field in the sample is given by: 



tH 



where 



J, 



H 



jk 



(11) 



(12) 



i.e. now we have a network with external Hall currents 
flowing into every node. Thus the paradigm of the 
hopping Hall conduction in the approach we adopt is a 
random resistor network with external Hall currents |l9| . 

To actually obtain explicit analytical results using the 
above procedure, one needs to know the distribution of 
Hall currents J/^ in the network. We will use ideas 
from the percolation theory picture of hopping trans- 
port 0, 0, ll^l to find this distribution. According to 
this picture, one can assume that all the hopping cur- 
rent in the network is carried by a percolating cluster 
of large conductances, such that for all conductances in 

the cluster G^J > Gc, where InGc {Tq/T)^/^ (we 

will limit our discussion to Mott's variable-range-hopping 
type conduction, but our theory, with minor modifica- 
tions, is applicable to other kinds of hopping transport 



<5'l23-F'l23 



G13G32 + G12G23 + G31G12 



(14) 



where (. . .)c denotes averaging over the on-site energies 
and intersite distances in the triad (1,2,3), such that the 
condition 

2r„/C + e„/r< (To/T)^/^ (15) 

is satisfied for every bond (ij) in the triad [l^ and 
5'i23 = sgn[h • (ri2 x ri3)] accounts for the orientation 
of the triangle with respect to h. Taking into account 
that the Hall current density is given by — (JxxE^ 
and / = axxL'^E, we obtain the following expression for 
the anomalous Hall conductivity 0,0|: 

S'l23-F'l23 



AH ^ T 2 

\ Lfl3(-f32 + tTl2Lf23 + tT3lUl2 



(16) 



Let us write down the averaging explicitly as an integral 
over intersite distances and on-site energies: 



1 



Oxy = Lalx^ I dri2dr23rfr3i / deidt2dfi 



X£'(ei)£'(e2)£'(e3) 



'S'l23-Pl23 



G13G32 + G12G23 + G31G12 



(17) 
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where is a normalization factor. The densities of 
states g(e) can be expanded around the Fermi energy 
g{e) ~ go + since hopping is restricted to a narrow 
interval of order r(To/T)^/^ around the Fermi energy. 
Substituting the expanded densities of states in Ea. l(T7|l . 
one finds that the zeroth order term, proportional to gQ, 
vanishes. The reason is that F123 changes sign under the 
inversion of all the on-site energies, while the bond con- 
ductances Gij are invariant under this transformation. 



Thus we obtain: 



AH 



La: 



2 

XX 



d In go / 5i23-Fi23(ei + £2 + £3) 



(18) 

where the average is now over a uniform distribution 
of the on-site energies. Approximating the average of 
the quantity inside the angular brackets by its maximum 
value and expanding to first order in h, one finally obtains 
the following expression for the anomalous Hall hopping 
conductivity: 



AH 
'xy 



Lat^4L_ ^To/Tf ,~^T.ITf'\ (19) 



?,/2 



where we have used the property that 1^/2 3> ii/2 in 
GaMnAs Analogous calculation in the case of the 

ordinary Hall effect gives 



OH 



' cet 



3/2 



(20) 



Thus, the ratio of the ordinary and anomalous Hall con- 
ductivities is 



AH 
tOH 



d In go 



de eB(- 



(To/T) 



1/4 



(21) 



Several concrete experimentally testable predictions fol- 
low from the above equations. First, according to p9|l . 

is proportional to the derivative of the density of 
states at the Fermi energy. This means that the anoma- 
lous Hall conductivity can be expected to change sign 
as the Fermi level crosses the density-of-states maximum 
in the impurity band. This is in contrast to the ordi- 
nary Hall conductivity cr^y^ which has the same sign ev- 
erywhere in the impurity band. Thus AHE can poten- 
tially be a very useful tool in elucidating the character 
of hole carriers in (Ga,Mn)As. The second prediction 
concerns the temperature dependence of cr^^ ■ As seen 
from Eq. H19|l , the leading temperature dependence of the 
anomalous Hall conductivity is oc exp[— (Tq/T)^/^], 
i.e. it depends on temperature in the same manner as the 
longitudinal conductivity (Txx- The same is true for the 
ordinary Hall conductivity. This, in particular, means 
that the Hall resistivity g^y ~ Qxx as a function of tem- 
perature, i.e. it diverges as T ^ 0. Thus our calculation 
is consistent with the view that the transverse resis- 
tivity of the "Hall insulator" diverges at T ^ 0. 



In conclusion we would like to point out that our the- 
ory seems to be consistent with recent magnetotransport 
experiments on digit ally- doped GaMnAs samples |2,1 |. 
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